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The formal expression for the most general polarization observable in elastic electromagnetic lepton hadron 
scattering at low energies is derived for the nonrelativistic regime. For the explicit evaluation the influence of 
Coulomb distortion on various polarization observables is calculated in a distorted wave Bom approximation. 
Besides the hyperfine interaction also the spin-orbit interactions of lepton and hadron are included. For like 
charges the Coulomb repulsion reduces strongly the size of polarization observables compared to the plane 
"l ] wave Born approximation whereas for opposite charges the Coulomb attraction leads to a substantial increase 

I . of these observables for hadron lab kinetic energies below about 20 keV. 
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I. INTRODUCTION 



Recently, Coulomb effects on polarization transfer from polarized electrons or positrons to initially unpolarized protons or 
1—1 ' antiprotons in elastic electromagnetic scattering have been studied in a distorted wave approximation at low energies fltl- These 
1-^, , studies were motivated by the idea to polarize hadrons by their scattering on polarized electrons or positrons in a storage ring . 
I ' However, in view of such a design it turned out that the considered observable, i.e. the total cross section for the scattering of 
Q initially unpolarized hadrons off polarized leptons to polarized final hadrons, the polarization transfer Pzooz, cannot contribute 
^ to a net polarization of the hadrons in a storage ring. The reason for that is that this polarization observable does not contain a 
Ph I genuine hadronic spin-flip process [11101, which is necessary for a net polarization change. Moreover, our previous numerical 
results were critizised by Milstein et al. ^ who had taken a partial wave expansion of the Coulomb scattering wave function 
instead of the integral representation used in ref. [1]. Indeed, it turned out that besides a minor error the main reason for the 
J> , gross overestimation of the polarisation transfer cross section was an accuracy problem in the numerical evaluation, namely, the 
' relevant quantity was calculated as a difference of two almost equal numbers multiplied by a huge factor 

For these reasons I have extended the previous study to the formal consideration of all possible polarization observables in 
2^ this scattering reaction including such spin-flip transitions using again the distorted wave Bom approximation. In addition to 
. the previously considered hyperfine interaction I have included also the spin-orbit interactions of lepton and hadron. In the next 
section the most general scattering cross section is introduced, defining the various polarization observables in terms of bilinear 
hermitean forms of the T-matrix elements. For the nonrelativistic form of the T-matrix with inclusion of hyperfine and spin-orbit 
interactions the detailed expression of the general scattering cross section is given, allowing for the polarization of all initial and 
final particles described by corresponding spin density matrices. In Section III I specialize to the case where the polarization 
of the final lepton is not measured to the so-called triple polarization cross section. For the numerical evaluation two different 
methods have been applied, a partial wave expansion as in ref. |4] and an integral representation of the Coulomb wave function 
according to ref. [7]. Results for the structure functions and spin-flip triple cross sections for the case of polarization along the 
incoming hadron momentum are presented in Section IV and a summary is given in Section V. Details for the evaluation of the 
hyperfine and spin-orbit interactions are given in an appendix. 



II. THE GENERAL DIFFERENTIAL CROSS SECTION INCLUDING POLARIZATION OF ALL PARTICLES 

Reviews on polarization phenomena may be found for lepton hadron scattering in [8], for nuclear physics in [9] and for 
nucleon nucleon scattering in [10]. I will consider hadron-lepton scattering in the cm. system, where hadron stands for proton 
or antiproton and lepton for electron or positron, 

h{p) + l{-p)^h{p') + l{~p'), (1) 

allowing for inital and final hadron and lepton polarization. The hadron initial and final three momenta are denoted by p and 
p', respectively. All possible observables of this reaction can be obtained from the "quadruple polarization" cross section for 
which the spin states of all initial and final particles are described by the corresponding general spin density matrices /^'/''(Pj^^ ), 
where the initial density matrices characterize the spin properties of target and beam and the final ones those of the detectors. It 
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is given by the general trace 



7 quadruple , 



O(P{,Pf,PJ,PI;0,, 



7^2^2(1 _,.|p/|)(l. 



ip;:i) 



Trace 



f^p\v{)p\v{)fp\Vl)p\V\) 



(2) 



where T = T{9, (p) denotes the T-matrix of the scattering process with {6, (j)) as scattering angles, and piV) the spin density 
matrix for a spin- 1/2 particle, with P characterizing the polarization of the corresponding particle in the initial and final states, 
respectively. The trace refers to the hadron and lepton spin degrees of freedom. The factor in front takes into account the final 
phase space, the incoming flux, and a normalization factor for the case of partially polarized final states. The invariant energy of 
the hadron-lepton system is denoted hy W — + Ei and the masses of hadron and lepton by and A//, respectively. In the 
cm. frame I use as reference system the z-axis along the incoming hadron momentum p. The x- and y-axes are chosen to form 
a right handed orthogonal system. 

In view of the fact, that in this work I am interested in the low energy regime, a nonrelativistic framework is adopted. The 
nonrelativistic density matrices for possible polarization of initial and final states of a spin-1/2 particle have the standard form 



p(P) = i(l + P-^) 



(3) 



with the vector P describing the polarization of the particle and a denoting the PauU spin vector One should note that in general 

From the basic equation ^ one obtains all possible polarization observables. In detail they are: 
(i) The unpolarized differential cross section: 



dao{9,(j)) 



0(0,0,0,0; 



(4) 



(ii) Beam, target and beam-target asymmetries of the differential cross section for unpolarized final states in the notation of 
Bystricky et al. ||i3l: 



dai: 



dn 



O(O,O,Pj,,Pl;0, 

dao{e,4)) 



dn 



1 + Pk,jAoojo {0, </>) + ^ Pl,kAoook (e» + E KA^oojk {9, < 

j k j,k 



(5) 



with the asymmetry vectors 



AoQOj{9,(l)) 



1 



d 



S° dP, 



-OiO,0,PlO;9,, 



da 



_pi ps 



2S° V dn dn 



25*" V dn 



f^o-p; ,-p; 

dh 



(6) 



(7) 



and the hadron-lepton asymmetry tensor 

.„ 1 



0(0,0, PL P?;0,< 
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(iii) Polarization of the final lepton or hadron for unpolarized beam and target: 

= j^^O{0,P{,0,0;e,(l)), (9) 

PjOOo{0,cl>) = J^^O(P{,0,0,0;^,</)). (10) 

(iv) Various correlations between the polarization of one outgoing particle and beam and/or target polarizations. For exam- 

ple, the outgoing hadron polarization for initial lepton polarization but unpolarized incoming hadron, the lepton-hadron 
polarization transfer is given by 

Pj00k{e,4') = ^— /^O(P{,O,O,Pj;0,0). (11) 

(v) Another interesting example is the hadron spin-flip of an initially polarized hadron by the scattering on an initially polarized 

lepton. It is a special case of the so-called "triple polarization" cross section with all particles polarized except for the final 
lepton as defined by 

=O(P{,O,P5,Pj;0,</>) (12) 



for the case P{ = ~PJi. i-e- 



O{-PlO,PlPl-0,ct>)- (13) 



This is the relevant quantity for the method of polarizing hadrons by electromagnetic scattering on polarized leptons in a 
storage ring L3n4J. 

A. The nonrelativistic T-matrix 

For the explicit evaluation of the trace in eq. (|2) one needs to know the spin dependence of the T-matrix. In a nonrelativistic 
approach but including contributions of the order AI~'^, the T-matrix contains the Coulomb, the lepton and hadron spin-orbit 
and the lepton-hadron hyperfine interactions. Separating the various contributions, the T-matrix is given in an obvious notation 
by 

f = fc + Tls, + Tls, + fss ■ (14) 

In detail, one has the Coulomb contribution 

Tc = 47raac , (15) 
the spin-orbit interactions of lepton and hadron, respectively, 

Tls,,, = 'inahi^f,-CT'/\ (16) 

and the hyperfine interaction 

fss = Anao-''- d -cr' , (17) 

•(->■ 

where d denotes a symmetric rank-two tensor, q = p ' — p the three-momentum transfer, and a denotes the Sommerfeld fine 

<-» 

structure constant. The tensor d can be decomposed into a scalar and a spherical tensor of rank two, i.e., a symmetric cartesian 
tensor with vanishing trace, 

d = + , (18) 
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where 

= do%, and do = -Trace(d) , (19) 
rf!? = '^u - doS^J ■ (20) 
Furthermore, the parameters ac, b;/;j, and d depend on what kind of approximation is used. These are: 

(i) Plane wave approximation (PW), corresponding to a pure one-photon exchange; the nonrelativistic reduction of the T-matrix 
including lowest order relativistic contribution reads 

Zh (2^11-1 2^ii\ Zi (2pLh~l 2^ih' 



8Mi V Ml 

WM/. „_ „ _ ^, (21) 



4:MlM, 

with P = p + p', and Z/j as the lepton and hadron charges, and /i; and /i/j as their magnetic moments, respectively. 
From this expression one reads off the parameters, keeping in the spin independent term the lowest order only, 

- (22) 

"i/h — ^<^i/h — — ' y^^) 

df^^ = c^«(g,?,-%), (24) 

where g denotes the unit vector along the three-momentum transfer q and q — |q |. The separation into a scalar and a 
traceless tensor according to (fT¥t reads 



2 
3 

,[2\PW _ SS 



^pw ^ _f_^ss^ (25) 



rfir - c^'m-M. (26) 



Furthermore, the strength parameters ai^e 



ep = ^f^^l_l + 2i^), (27) 



4Mi V Ml M,_ 

One should note that the strength parameter of the hadronic spin-orbit interaction is about three orders of magnitude 
smaller than the parameter of the leptonic one, because their ratio is approximately 

4J«2/x;,^«3-10-3. (30) 
cf^ Mh 

(ii) Distorted wave approximation (DW) using nonrelativistic Coulomb scattering wave functions 



sinh 7r?7c 

and 



V'p ^+^1- ) = e-*"^ e'P-- lF,{-^r,c, 1; ^{pr - p • r )) (31) 



<(-)(r) = (^^W(r))*, (32) 
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where i^p^^^ denotes the incoming and outgoing scattering waves ifTlll . respectively. Here, iFi(a, denotes the 
confluent hypergeometric function. In the expression for the scattering wave in eq. dSTT l I have akeady separated the 
constant Coulomb phase factor e"'^ with 



ac = arg[r(l + irjc)] 



(33) 



because it will disappear in the observables. The relevant quantity for Coulomb effects is the Sommerfeld Coulomb 
parameter 



rjc = aZiZh/v 



with V denoting the relative hadron-lepton velocity. 
Within this approach one finds 



u,DW 



e^^^ag^ , with 0c (e) = -?7cln[sin2(0/2)] , 



r.LS 



d. 



DW 



47r / r'i P 



Separating again the hyperfine contribution into a scalar and a traceless tensor, one obtains 



One should note that the tensor df^ in eq. (l37T i is symmetric as well as d'^' 



(34) 

(35) 
(36) 

(37) 

(38) 
(39) 



B. The general scattering cross section and polarization observables 

Evaluation of the trace in eq. (|2]l yields the following general expression 



O(P{,P/,Pj,P?;0,. 



a,[3£{C,LSi,LSh,SS} 



where the various contributions are defined by 



MfMl 



^2M/2(l+|p/|)(l+|p/|) 



Trace 



with Tq, defined in eqs. (fT5l)-(fT7l). One should note the relation 



(40) 



(41) 



(42) 



from which follows that Sa ■= Sa,a is real. 

Separating the diagonal contributions (Sa) from the interference terms {Sa,i3 for a ^ (5), one obtains for the "quadruple 
polarization" cross section 



dn 



JG{C,LSi,iSh,SS} 



a<l3e{C,LSi,LSH,SS} 



(43) 
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Explicitly, one finds in terms of the different contributions to the T-matrix in eq. ( fT4] i for the diagonal terms 



Sc{0,' 



Vo\ac\^n+U+ , 
Vo n+ {U- br -hi +2 Re [{hi . Pf ) * (b, . Pj)] ) 



^o(n,;nri?o-P,;-G-P 

Im 

jst 



(44) 
(45) 
(46) 



+2 Re 



[U-P{ -D-Pl + il^h) 
(n- {P- • H) + 2 5] Pi:X.sPb E,st) h) 

jst 

{p{ ■ T ■ Pi) {p]-d- PD + (Pi • T ■ Pi) {PI • d • pf )]) 



where 



Vn 



W^{l + \Pl\){l 
The following quantities depend on the polarization parameters 



|P{| 



p± 

Qh/i 



i±p{„ -P 



h/l 



h/l ^ h/l I 

- iP{„ X P 



P/i/i ^ ^h/i 



h/l 



h/l 



(47) 



(48) 



(49) 
(50) 
(51) 



Furthermore, I have introduced for convenience the following quantities, which depend on the hyperfine interaction tensor , 



Do = 

Ejst = 

Gij = 

H, = 



fc 

, ^ 1 

2_^d*jdj, = -Trace(£)), 

ij 

kl 

^ ^ ^ils ^jrat dij^dst ^ 
Imst 



(52) 
(53) 
(54) 
(55) 
(56) 



klm 



where Siti denotes the totally antisymmetric Levi-Civita tensor in three dimensions. These quantities are functions which depend 
on the scattering angles {6, </>). 

Correspondingly, one finds for the interference terms 



Sc,LSi ( 

Sc,ss{ 
Sls,,ss{ 

SLSh.Ssi 



= Voa*cIl+hi-Qi, 

— Sc.LSl {S, (t^)\h4r^l , 

= Vo a}j Qh ■ d ■ Qi , 

= Vo {hi ■ Qi)* {hh ■ Qh) : 



voQh-d- in- K - *br X p- + (br • PDPf + (br • pf)p: 



s 



LSi,SS 



{0,^)\h 



^l 



(57) 
(58) 

(59) 
(60) 

(61) 
(62) 
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Using the separation of the tensor d into a scalar and a traceless symmetric tensor according to eq. JlSt , one finds 

Do = 3|doP+E4f41. (63) 

i,k 

= 2|do|'%+4i?e(d*dg') + A(f , (64) 
= 2i?e;^4fdg, (65) 



where 



Ejst 



Mop + E ^i^^i 41 + ( E ^^-^^ ^s4' - ^ 0*) , (66) 

kl I 

G„- = 2|doP%-2i?e(dS4?) + Glf , (67) 

i/. = Y.^^Hdf:^d^Z, (68) 



Gg) = ^ . (69) 



Imst 



It is now easy to see that the vector H is purely imaginary and that the tensor Gij is real and symmetric. Furthermore, one notes 
the symmetry property 

E*,, = ~E,ts ■ (70) 

It suffices to evaluate the spin-orbit vector b and the hyperfine tensor d for = 0, because then the values for an arbitrary can 
be generated by a rotation around the z-axis exploiting their rotation properties. Examples are given in the following section. 



in. THE TRIPLE POLARIZATION CROSS SECTION 



I will now specialize to the case where only the final lepton polarization is not analyzed, i.e. Pf = 0, but all other particles 

are completely polarized (|P^^''^| = 1, |P;| = 1). This case is of particular interest for the polarization transfer in a storage 
ring L3,^. The corresponding "triple polarization" cross section has the form 



= o(p{,o,p5,p: 



E 

a<l3e{C,LSi,LSh,SS} 



In this case, the lepton polarization quantities in (|49] l through (ISTT i become 



and one finds for the diagonal terms 

\9) 



^triple f 



gtri'ple / 
LSi \ 
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nf = 1, 

VTVl\acm\ 

V (n- b;, • hn 



Q; 



1 1 
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P. • H + n, (PI • H) + 2 ^ KAA Ejst 



Im 



(71) 

(72) 

(73) 
(74) 

(75) 
(76) 



jst 



and for the interference terms 










ill) 




= Vachh-Q,h-, 


(78) 




= Va*cCth- d-P], 


(79) 
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— V (Yi, ■ P?")* Chi, • Ou\ 
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= VPld-(ll^K-ihlxP^ + (hl- Fl) Pi + {hi ■ Fi) Pi) . 


(82) 


where 










(83) 



From now on as a further specialization, I will consider only polarization along the incoming direction which is chosen as 2-axis. 
Then with 



= Af5, P± = (Al±A^)5, and n± = 1 ± A^A^ , 



(84) 



where A^ = A^j ± A^, one obtains 



dO. 



= (1 + a5,a{) [sc{e) + Soie, 0) + L\,{e, 0)] + (i - a5,a{) <f) + K {si{e, <f>) + Ll{e 



+a: 



's,{e, 4>) + yi(s^{e, <f) + 4^(0, </.))] + a+ [a^ (52+(0, </.) + 0)) + L\{e, 4> 



+AlA{[s2(^,(/.) + Ai53(^,(/>) 
The diagonal contributions are 

Sc{e) 

and the interference terms 



V\ac\\ 

v(3\dof + J2Dni0,<l^)), 
3 

iV^e3jkDjkiO,cl)) = i{Di2{e,ct>) - D2iie,cP)) = -2 7m(£»i2(^,<^)) , 

jk 

2v(^Re(d*Ji{e, <t>) ~ dfhe, (i>rSi{e, 4,)) - \do\' + \dfl{0, </>)p) , 

2v{2Re{d*od§{0,<P)) - 2|do|2 - - i^22(^,</')) , 

-4vim(d^il{e,ct>rSil{e,^)), 



{0, 4>) &r/,,2(^, ^) - 41 {0^ <t>) , 



= 2VIm 



S+ {0, 4>) = 2V Re a*c{do + 41 {0, <^)) 



(85) 

(86) 
(87) 
(88) 

(89) 

(90) 
(91) 
(92) 

(93) 
(94) 
(95) 



Here I have already used the fact that the spin-orbit vector hi/h is perpendicular to the ^;-axis, i.e. its third component vanishes. 
Indeed, as shown explicitly in the appendix, it has the form 



hi/h{0,^)=ib'J\0){-sm^,coscP,O)=ib'o^\0) 



P X p 
|p' X p| 



(96) 
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However, not all of the contributions in eqs. (l8&t through (|95l l are nonzero. In fact, I now will show that the diagonal 
contributions Si and S3 and the interference terms l'/'' vanish identically. To this end I first will consider the (/)-dependence of 
the tensors cfi^^{9, (f>) and D{9, 4>). It suffices to evaluate them for </> = from which one obtains S'^\9, (p) and DlO, <j)) by a 
rotation around the z-axis by an angle Introducing the notation d'^l'^ (6') = d''^) (6*, — 0) and correspondingly I?° = D{6,0), 
and using furthermore the fact dfj — ~ ^'^^ finds 



cos^ (t> + 4f °W sin' 5(4t°W - 4f °W) sin 20 4f °(^) cos A 



i(4t°(e) - 4f sin 20 4t°W sin' + 4f cos' 4f °(^) sin0 



V 



4f°(6') COS0 



4f "(f) sin0 



j[2],0 
^33 



(97) 



{0) / 



From this representation one notes immediately that 

,\2 



4i(^.<^)*42(^>'^) = |rfit"W|cos0sin0 



is real and thus S3 



( 



D°{9) 



i(^,'^) = l4f° 
according to eq. d92] i. Furthermore, with 

i4t"wp + i4f °(^)i' 4t°w*4f + 4f °(^)*41'°w\ 



(98) 



Krm 

V4f "(^)*4t°w + 4f °w*4f "(^) i4f "wp + i4f 

which formally has the same structure as dP^'^{6), one finds 

(0)cos2 + ^02(0) sin' i(i??i(0)-i?O2(f?))sin20 7^03(0) cos A 



(99) 



Here again one notes that 



l{D°^{e) ~ D^^i^)) sm2(t) D°^{9)sm^(j) + D^2i&)co8^(j) D°^{e) sin (j) 
V D%{e) cos ct> D%{e)sm4, D%{e) J 



Di2{e,cl,)^-{DUe)-DU0))sin2<j, 



(100) 



(101) 



is real because according to eq. (|99T l [9) and Z?22 (^) '"^^l ^'^'1' therefore. Si = 0. For the spin-orbit interaction vector 
bw/j one finds from eq. ( |96] l |b;//j(6', </>)p = |foj/''(6')p, which means that ij/'* is ^-independent. Moreover, using eqs. ( |96] l and 
(|97]l one obtains 



(6*, (j>) = 4t/ i?e(dfj'°(6') cos (t){~b'J'\9) sin 0) + sin ^^/''(e') cos 



0. 



(102) 



Finally, with the help of eqs. (|97T l through dlOOt one finds that the remaining structure functions become (^-independent, which 
is easy to understand since all polarizations are assumed to be along the z-axis ruling out any 0-dependence. Thus also the cross 
section simplifies to 



= (1 + Aj^A^ \Sc{9) + So{9) + 4(0)1 + (1 - XlXi)L^oiO) 



+A5 



XJ;(S^{9,^) + L'^{9)) + X+(s+i9) + L^,{9))] + ^l^iS2{9) 



(103) 
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where the structure functions are given by 



L. 



i/h. 



So{0 



5^ 



= 2VRe 



V 



(3Mo| 



2VRe 
2V 



33 



j[2],0 

hi 



2|4f°(0)P 



(0)*4f°(0) 



Mo I 



2V 



2 Re 



{d*dfl\e))-2\do\'-\dri0)\' 



l4f°(^)|^-Mg'°WI^ 



(104) 
(105) 

(106) 

(107) 
(108) 
(109) 



The expression in eq. (1103b is an extension of the triple polarization cross section given in eq. (16) of ref. ^ by including the 
contributions from the hadron and lepton spin-orbit interactions. (As a sideremark: The corresponding cross section for final 
lepton polarization is obtained from (l85t by the substitutions A{ X( and AJ^ -o- A|.) 

As a special case, one considers the so-called hadronic spin-flip cross section for complete hadron polarization, i.e. A^^ = 
— A{ = A/i = ±1, and the non-spin-flip cross section with AJ^ = A^ = A/, = ±1. For the spin-flip cross section one finds 



2L^{9) - S2{e) + 2XIX,, S^{9, 4>) + L'^{9) 



(110) 



It is governed by the hyperfine terms 52 and 5*2 and the hadronic spin-orbit interaction via Lq and On the other hand, the 
non-spin-flip cross section is given by 



do-"*** 



dn 



= 2 



Sc{0) + So{0) + L'oie) +S2{0) + 2XlXh S+{9,(b)+Li{e) 



(111) 



Its helicity independent part is overwhelmingly dominated by the Coulomb term Sc with additional tiny contributions from the 
hyperfine and the leptonic spin-orbit interactions. The difference of the non-spin-flip cross section for Xh = ±1 



2 V dn 



dn 



= 2A 



(112) 



has been considered as lepton-hadron polarization transfer in This polarization transfer is dominated by the hyperfine 
structure function S2 because the additional spin-orbit contribution is comparably small as shown in the next section. It 
differs by a factor 2 and the presence of L2 from the polarization transfer Pzqoz for the scattering of unpolarized hadrons on 
polarized leptons as considered in ref. [ 1 ], where I had considered only the leading term S*^, the interference between Coulomb 
and hyperfine amplitudes, neglecting higher order contributions. The more complete expression reads 



daoi 



da 



zOOz - 



triple ^ 



dn 



dxidx] 



dn 



-L'' 



(113) 



In addition to 5^, it includes S2, which is quadratic in the hyperfine amplitude, and L2 and L2, the contributions from the 
interference of hyperfine and leptonic and hadronic spin-orbit amplitudes, respectively. However, the largest of these additional 
terms, L2, is still quite small if not negligible compared to 5^. 



IV. RESULTS FOR STRUCTURE FUNCTIONS AND POLARIZATION CROSS SECTIONS 

For the evaluation of the structure functions in eq. ( 1103) and the corresponding cross section I have used both methods for the 
calculation of Coulomb distortion, the integral representation as well as the partial wave expansion. The integral representation 
has been used mainly in order to check the convergence of the partial wave expansion as is described in detail in the appendix. 
Thus all results presented in this section are based on the partial wave expansion (PWE). For the hyperfine amplitude it was 
found that an expansion up to a partial wave with Imax — 2000 was sufficient, but for the spin-orbit interaction, beeing much 
slower convergent, Imax — 4000 was taken. 
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FIG. 1: Structure function So{6) of the triple polarization differential cross section in plane wave (PW) and distorted wave approximation for 
lilie charges (DW^) and opposite charges (DW~) for various proton lab kinetic energies Th . 



A. The structure functions 



First, 1 will discuss the various structure functions which determine the triple polarization cross section of eq. (11031) . For the 
plane wave approximation, one finds easily the following expressions 



L 



rl/h, 

^0 I 

^o( 



S+{0 

S2{0 



v{4/%rcot'{e/2), 



LS „SS 



2V{c' 
VcSS 



cos2(0/2), 



SS\2 



cot2(6'/2) 



2p2 

= -2F(c^^)2(l + sin2(0/2)) 



(114) 

(115) 
(116) 

(117) 

(118) 
(119) 



The structure functions, evaluated in the cm. frame for several lab kinetic energies, are shown in Figs. [T]through|6]for the various 
approximations, i.e. plane wave approximation (PW) and Coulomb distortion for like (DW+) and opposite charges (DW^). 

The diagonal structure function 5o in Fig.[T] induced by the hyperfine interaction, shows a rather flat, almost constant angular 
behaviour. Its size scales roughly proportional to the inverse of the kinetic energy Th- Compared to the plane wave approximation 
(PW), the distorted wave approximation is strongly enhanced for unlike charges (DW~) and strongly suppressed for like charges 
(DW+) by several orders of magnitude for the lowest lab kinetic energy of Th = 0.001 MeV corrsponding to rjc ~ 5. This 
enhancement resp. suppression is increasingly reduced with growing kinetic energy T/, and approaches the plane wave result 
above Th — 10 MeV. The pure Coulomb contribution Sc, however, is much larger by more than ten orders of magnitude. 

With respect to the other two diagonal structure functions from the leptonic and hadronic spin-orbit interactions Lq and Lq, 
it suffices to show only the former one in Fig. |2] because Lq differs in magnitude only by the factor {c\sl ^^ls)"^ being smaller 
by about five orders of magnitude. One readily notes that Lq exhibits a strong peaking in the forward direction only and tends 
to oscillate at small angles for the lowest Th considered here. Over the whole angular range, especially in the forward direction, 
Lq is much larger by several orders of magnitude than 5o but it is still almost negligible compared to the size of Sc- The effect 
of Coulomb distortion is qualitatively similar to what one observes in Sq. 
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Only these diagonal structure functions contribute to the unpolarized cross section. However, as already mentioned, their 
relative contribution is extremely small as can be seen by comparison with the pure Coulomb structure function Sc which is 
also shown in Fig. |2]indicated by the large reduction factor applied to Sc- 

The two hyperfine-hyperfine interference structure functions S2 and S2 , shown in Figs. [3] and 21 which both are negative 
througout, exhibit a similar pattern, a smooth angular distribution with a slight decrease in size at backward angles for the two 
lowest kinetic energies, but with a slight increase for the two higher kinetic energies like the PW result for all four cases. Again 
one notes sizeable enhancements for opposite charges by Coulomb distortion and suppression for like charges for the lowest 
energies Th. The distortion effect decreases with increasing Th- Also these two structure functions are quite small like Sq 
because they are quadratic in the hyperfine amplitudes. 

Much larger is the third interference structure function S*^ because it is an interference between the Coulomb and the hyperfine 
amplitudes. Thus it is strongly forward peaked. For this reason, it is displayed in Fig. |5] multiplied by sm^{9/2). Moreover, 
Coulomb distortion induces a strong oscillatory behavior, again in conjunction with a large enhancement for opposite charges 
and strong suppression for like charges. As expected, the distortion effect diminishes with increasing kinetic energy Th- 1 would 
like to point out that in the corresponding Figs. 1 and 2 of ref. yj] erroneously a factor sin^(2^?) instead of the indicated factor 
sin^(^?/2) has been applied. 

Finally, the spin-orbit-hyperfine interference structure function Fig.|6]is comparable in size to 5*2 and 5*2^ but exhibits 
quite a different pattern. At the lowest kinetic energy Th it is strongly enhanced by distortion for opposite charges and possesses 
a pronounced broad minimum around 100°. It falls off at forward and backward angles with many oscillations in the forward di- 
rection. With increasing kinetic energy the minimum moves towards smaller angles with fewer oscillations. Again the distortion 
effect decreases strongly with increasing energy. Like the diagonal structure function Xq, the hadronic interference structure 
function is quite negligible. 



B. The triple polarization cross section 

Now 1 will discuss the triple polarization cross section of eq. ( |103t . Previously, in ref. [3] only the hyperfine amplitude besides 
the Coulomb one has been considered, whereas the hadron spin-orbit interaction has been included akeady in ref. [4]. However, 
as mentioned above, its contribution to the helicity dependent part of the spin-flip cross section in eq. (II 10b is negligible, whereas 
in the helicity independent part the diagonal contribution Lq is comparable in size to 5*2 in the forward direction (compare Figs. |2] 
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FIG. 7: Absolute value of spin-flip cross section da^ / dO, for initial hadron polarization parallel to lepton polarization along the initial relative 
momentum in plane wave (PW) and distorted wave approximation for like charges (DW"*") and opposite charges {DW~) for various proton lab 
kinetic energies Th. 
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FIG. 8: Absolute value of spin-flip cross section dcr^J /dQ for initial hadron polarization opposite to lepton polarization along the initial 
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30 60 90 120 150 180 30 60 90 120 150 180 




Tj, = 0.01 [MeV] 



30 60 90 120 150 180 30 60 90 120 150 180 

e [deg] 
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momentum with the hyperfine interaction alone (DW(hfs)) and with inclusion the of spin-orbit contribution (DW) for like charges (DW"*", left 
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andO. Much more important is the leptonic spin-orbit contribution which, however, appears in the non-spin-flip cross section 
only (see eq. dl 1 lb ) where it is buried completely by the Coulomb contribution Sq- 

The results for the spin-flip cross section for parallel initial spin orientations of hadron and lepton is shown in Fig. |7] while 
the one for the opposite spin orientation in Fig.|8] One notes again the strong influence of Coulomb distortion. Furthermore, the 
leptonic spin-orbit interaction plays a relatively important role in the region of the minimum as can be seen in Figs. |9] and [TOl 
where a comparison is exhibited with the case for which the spin-orbit interaction is switched off (curves labeled "(hfs)"). One 
readily notes a substanial increase when the spin-orbit part is included compared to the pure hyperfine case. Furthermore, at the 
lowest two energies the spin-orbit interaction induces oscillations, in particular in the forward direction. The difference of the 
two spin-flip cross sections determines the net hadron polarization in a storage ring of initially unpolarized hadrons scattered at 
polarized leptons. 

The non-spin-flip cross section in eq. (II 1 lb is overwhelmingly dominated by the Coulomb contribution Sc- The small 
dependence on and leads to different scattering strengths for the hadron polarization parallel or antiparallel to the lepton 
polarization (see eq. (II 12b ). The resulting lepton-hadron polarization transfer PzOOz is shown in Fig.[TT] The dominance of the 
hyperfine amplitude is clearly seen. 



C. The integrated structure functions and cross sections 

Finally, 1 will present results for the integrated structure functions and spin-flip cross sections which are the relevant quantities 
for the polarization build-up in a storage ring. They are defined by the integration over the solid angle except for the small 
cone in the forward direction with < 9,nin, where the minimal scattering angle is defined by the requirement that the impact 
pai^ameter should not exceed a given value b 

9,nin = 2 arctan(7yc'/0 i (120) 

with I = bp as classical angular momentum. In the present work 1 have chosen b = 10^'^. The choice of this value has been 
justified in ref. tSJ. The dependence on this parameter is discussed below. Thus for any structure function or cross section O{0) 
I define 



(O) = 27r / d{cos0)O{e). (121) 
For the plane wave approximation, one finds easily the following expressions 

{l'J'\0)) = -27r\^(cf/f)2(ln(sin(0™„/2)) + icos2(0„„/2)) , (122) 

= -27rFcf/fc^^(l-sin2(0„,„/2)(l + cos2(0„,„/2))) , (123) 

{So{0)) = 87TV{c^^f, (124) 

{S+{0)) = 4^F— [ln(sin(0™„/2))--cos2(^^„„„/2)j, (125) 

{S^{0)) = -47rF(c^^)2cos2(0,„,„/2)(l + sin2(^?„„„/2)), (126) 
{S2i0)) = -47rF(c^^)2(3-sin4(0„„„/2)). (127) 

Thus in plane wave all of the integrated structure functions except for (5^) are almost independent of T/, except for a very weak 

dependence via the minimal scattering angle. One should note the logarithmic divergence for 0min — > in and (5^). 

It corresponds to the logarithmic divergence in the angular momentum / of the partial wave expansion noted in ref. [4] which 
appears when integrating over the whole range of scattering angles. 

The results are exhibited in Figs. [T2]through[T4lfor the structure functions and in Fig.[T5]for the spin-flip cross sections. The 
integrated structure functions show a strong increasing influence of Coulomb distortion with decreasing hadron kinetic energy 
T/j leading to large enhancements for opposite charges and strong suppression for like charges compared to the plane wave case. 
The only exception is (S*^) for which Coulomb distortion results in a reduction for both cases, however much stronger for like 
charges. The reason for this feature is the strong oscillatory behavior in the angular distribution of (5^) at low Ti^ (see Fig.|5])- 

The integrated spin-flip cross section is given by 



(ai^)=2(iJ>-(52)+2±[(52-> + (4) 



(128) 
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The corresponding integrated spin flip cross section of Milstein et al. |4] reads according to their eq. (21) 

«0M55 = ^(^)'[(27rr;c)'(y-ln2) + ln(/„,./,7c)2T(27r,7c)'] ■ (129) 

The appearing logarithmic divergence in the angular momentum I is regularized by choosing a finite Imax determined by the 
classical relation Imax = bp which corresponds to the choice of a minimum scattering angle in the present work. The integrated 
strength of the spin-flip cross sections in Fig.[T5]shows as expected with decreasing Th a growing strong influence of Coulomb 
effects via the hyperfine and hadronic spin-orbit interactions. The latter is only important in the spin-independent part of the 
spin-flip cross section while its influence in the spin-dependent part is negligible. The results of Milstein et al. (01, shown in the 
upper panels of Fig. [T5]for opposite charges, are comparable to our results but display a slight overestimation which is probably 
due to different approximations in [4]. The dependence of the integrated spin-flip cross section for opposite charges on the 
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100 



regularization parameter h is exhibited in Fig.[T6]for h — 10^, 10^° and 10^^ fm. It appears to be quite weak. 

The relevant quantity for a polarization build-up in a storage ring is the ratio of the spin-independent part over the spin- 
dependent part 



{3-2 



(130) 



which is shown in Fig. [T7] . One readily notes the reduction of this ratio by the hadronic spin-orbit interaction, in particular 
quite strong at higher energies but only about 12 % at the lowest energy. This fact clearly shows the importance of the hadronic 
spin-orbit interaction besides the hyperfine contribution. 

Finally, I show for completeness in Fig.[T8]the integrated polarization transfer cross section (PzoozCo)- It is dominated by the 
hyperfine interaction whereas the spin-orbit contribution is negligible. 
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FIG. 17: Ratio of the spin dependent part of the spin-flip cross section over its spin-independent part as function of the proton lab kinetic 
energy Th for plane wave approximation (PW) and with Coulomb distortion for opposite charges (DW^): present calculation and the result of 
ref. 0] (MSS). For the curves labeled (hfs) only the hyperfine amplitude is included. 
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V. CONCLUSIONS 



Formal expressions for polarization observables in electromagnetic hadron-lepton scattering have been derived within a non- 
relativistic framework including the central Coulomb force as well as the lepton and hadron spin-orbit and hyperfine interactions. 
The latter have been treated in a distorted wave approximation. Special emphasis has been laid on the triple polarization cross 
section with polarizations of the initial hadron and lepton and of the final hadron along the incoming hadron momentum. The 
structure functions which determine the differential triple polarization cross section have been evaluated in plane and distorted 
wave approximations for hadron lab kinetic energies between 1 keV and 100 MeV. 

For the evaluation of the Coulomb distortion two different methods have been employed: (i) an integral representation of the 
nonrelativistic Coulomb scattering wave function and (ii) a partial wave expansion. These two independent methods have served 
as a mutual check for the numerical accuracy of the results. 

As expected, the distortion effects are very important at low energies in the small polarization observables, which are driven 
by the spin-orbit and hyperfine interactions, leading to sizeable enhancements for opposite charges and suppressions for like 
charges according to the Coulomb attraction resp. repulsion. This is shown in detail for the structure functions of the triple 
polarization cross section and for the special case of the spin-flip differential cross section. 

The leptonic spin-orbit interactions plays a minor role in the non-spin-flip cross section in its spin-dependent part, which, 
however, as a whole is by many orders of magnitude smaller than the spin-independent pait, dominated by the Coulomb term 
Sc- The influence of the spin-orbit and hyperfine interactions on the unpolarized cross section is almost negligible in the whole 
range of energies studied here. 

With respect to the integrated spin-flip cross sections our previous work has been extended by the inclusion of the hadronic 
spin-orbit interaction which shows a non-negligible effect in the spin-independent part changing sizeably the ratio of the inte- 
grated strength of its spin-dependent over the spin-independent part. 
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Appendix: Evaluation of hyperfine and spin-orbit interaction 

For the evaluation of the Coulomb distortion of the amplitude b of the spin-orbit interaction in DWBA as given in eq. (|36] | 
and the tensor amplitudes d^^I of the hyperfine interaction in eq. ( |39] l two different method have been applied: (i) An integral 
representation and (ii) a partial wave expansion of the Coulomb wave function. For convenience I set in this appendix rj — rjc- 
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1. Integral representation 



In iH] a detailed description of this method for the evaluation of d'^l has been given. Therefore, I will only summarize the 
result. The method is based on an integral representation of the confluent hypergeometric function as proposed in ||7| 



with 



Q(,).^,and/(t,,) = e--A 



(A.l) 



(A.2) 



[21 

With the help of this representation, the hyperfine tensor d\j and the spin orbit vector b are expressed as two-dimensional 
integrals, as described below. 



For the hyperfine tensor one finds 



Hyperfine interaction 



where 



TTT] 



is a normalization factor and 



dlfit,ri,0) - A,,{t,l,r^,0)lssic{t,l)) 
^ dt' 



Here I have introduced the tensor 



with 



and 



Furthermore, Iss (c) denotes the integral 



i ^e-''>"(i-*')(A,,(t,l)/(c(t,l))-e''""*'A,,(t,0/ss(c(t,i'))) 

Aij{t,t') = 3di{t,t') dj{t,t') - Sij , 
-a +'\ Pt-p't' 

g{t,t') = [i2+t'2_2tt'cOS0]l/2. 

1"°° dx 

Iss{c) = I —e^'^hix) 



1 1 2 1 

3^2^ ~ 4 



with Q{x) as the Heaviside step function and 



c{t,t') 





c+ 1 






c- 1 


-t- 


t' 



i7re(l -c)) , 



9{t,t') 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



(A.8) 



(A.9) 



(A. 10) 
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One should note that d\j and d^j (t) are functions in 6 and 0, the scattering angles in the cm. frame. However, as mentioned 
above, it suffices to choose = 0. 
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The remaining integrations over t and t' in eqs. ( IA.3l l and (IA.5l l are evaluated numerically. Details are presented in |[T|]. 
However, I would like to mention that in contrast to eq. (A. 36) of ref. 1 1], I found it more advantagous to use for the numerical 
evaluation a transformation of the integration variable y = — ln(l — t) resulting in 



^ dt 



1 - t 



^^'-('-'\g{l)-git)) = / dye^^y{gil)-g{l~e-y)). (A.ll) 



As mentioned in the appendix of ref. [jj, it is useful to separate the contributions of the real and imaginary part of the radial 
integral Iss for 77 > according to 

df;{rj,e) = 7V(77)(7ef/(7y,0)+Zf (r;,0)), (A.12) 

where TZ/X^^ - from now on called reduced amplitudes - refer to the contributions of Reiss and Imiss, respectively (see 
eq. ( |A.9t ). to djf' (t, -q, 9) in eq. ( |A.5t . With respect to the symmetry under sign change of rj, one finds easily the properties 



from which follows 

7[2] 



7^f/ (0,-77) = (7^f/(77,^?))^ (a.13) 

Iff{9,-v) = -{iffMY, (A. 14) 



,-77) = iV(-r;)(7^f/ (77,0) -1^-^(77,0)) , (A.15) 



Thus it suffices to evaluate the reduced amplitudes TZ/X^^ for 77 > from which one can determine (i[^'(0, ±7;) applying 
eqs. (IA.12b and ( lA.lSb . On the other hand, given (i^j {9, ±77) the reduced amplitudes are obtained from 



l<-^3^^^^) 2V 7V(77) + iV(-77) ) 2N{t^)\ '^^^' '^ ' ^' )' 

(^''^^ - 2\-lm m:^)-2m)y ''^^' ^~ d^A^^^-^) )■ (A. 17) 



(A. 16) 



These relations are useful for the comparison with the partial wave approach. In addition, they show that for large positive 77 the 
second terms in eqs. ( |A.16t and ( |A.17t are strongly suppressed by the factor e^^'"'. Consequently, one finds that TZfj^{r], 9) « 

lfj^{ri,9) for positive 77 S> 1. This feature poses a serious problem for the numerical evaluation of df^{9,~ri) for 7; ^ 1 
because it is proportional to the difference of the reduced amplitudes and multiplied with a huge number (iV(— 77)) according to 
eq. dA.lSt . Thus with increasing absolute value of 77 more and more significant digits are lost in the difference for negative 77. 
An example is presented later. 

a. Spin-orbit interaction 
Following the analogous steps for the spin-orbit interaction, one finds 

b(0,0) = 7 6o(7?,e)^7^ (A.18) 
IP X p I 



where 



60(77,0) = c^^sin07V(77) [60(0, 1,^7) r^e-"'i"(i-*)(6o(e,l,^7)-e"'^"46o(0,t,^)), (A.19) 



- g{tAY +'"^1 l-f' V .g(i,l)2 git,t'r ^ 

with g{t, t') in eq. ( IA.8I ) and the radial integral 



(A.20) 



dx 



c 







c+1 



Ils{c) ^ — e*^^ji(x) = l--ln ^ +7— e(l-c). (A.21) 



c- 1 



, TTC 



2 
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The numerical evaluation is analogous to the one for the hyperfine interaction. For 77 = one finds bo{0;0) = 
c^'^ sin6'/(4sin2 6/2) and thus 

b^^(0, 0) = loot (0/2) c^^^H;^ (A.22) 
/, IP X P I 

in agreement with eq. (l23T l. 

Also in this case I introduce reduced amplitudes by separating the contributions of the real and imaginary parts of the radial 
integral I^s acording to 

boiv,9) = iV(r,)(7^^^(r;,0) +1^^(77,0)) , (A.23) 

such that 

bo{e,-r,) = iV(-r,)(7^^^(77,0) -1^^(77,0))*, (A.24) 

and 

^"^'{^,0) - ^j^{bo{v,0) + e-'^no{e,~rjr), (A.25) 
I'-'iv^O) = ^^{bo{v,e)~e-^-no{9,-rjr). (A.26) 
Again numerical problems arise for negative 77 with I77I S> 1 as for the hyperfine interaction outlined above. 

2. Partial wave expansion 

The expansion of the Coulomb wave function into partial waves reads 

A-TT 

V'^+Hr ) = — E Fiirj,pr)Y;^if) YUp) (A.27) 

^ l,m 

where the radial function F; is given in terms of the confluent hypergeometric function iFi (a, 5, z) 

Fi{r^,p) = Ci{jj)e'P p'+\Fi{l + l + i7^,2l + 2,-2ip), (A.28) 
and constants depending on the Sommerfeld parameter 77 

Coin) ^.2 



(2; + l)!! 



Diiv'), (A.29) 



Co(77) = e-t"|r(l + z77)|=e't'^ /_Z!L^, (a.30) 

sinn(7r77j 



In the above expression, I have separated as in eq. (ISTT l the I ~ phase cto = ac for convenience, coinciding with the Coulomb 
phase ac given in eq. ( |33] |. The remaining partial wave phase is given by 

ai=ai-ao, (A.32) 

where for I > 



25 



Evaluation of the various contributions to the scattering matrix as listed in eqs. ( |36] | and ( [39] l leads to the following expressions 
which are still operators in spin space 



1=1 



ij l'=0 1=0 

where I have introduced for convenience in the notation of Fano and Racah lfT2ll for irreducible spherical tensors 



[K] 



1 X 



(A.34) 
(A.35) 

(A.36) 
(A.37) 



The coefficients are given in terms of the radial matrix elements Rin (note the meaning of the"hat symbol": / — \JH + 1) 



/' I 2 




Ri'. 



with 



I'l 



dr 



Jo J' 



Fi'{r],pr)Fi{r],pr) = 4 



P- 



dp , . , . 
— Fi'{r],p)Fi{'q,p) 



(A.38) 
(A.39) 

(A.40) 



Besides the radial integral Ru, only Ri'i — Rw for \l — l'\ — 2 are needed in view of the selection rule of the 3j-symbol in 
eq. ( IA.39b . This radial integral is well known in Coulomb excitation (see e.g. fi3\). The explicit form of the radial integral is 
also given in |i4J. For = / and I > one has 



Ru — 



1(1 + 1) 



fiiv) 

21 + 1 



with 



fiiv) 



sinh (ttt]) 

One should note that /; vanishes for t] — 0. For \l' — l\ = 2 one has 

RlI+2 



k=l 



3\l + 1 + iri\\l + 2 + iri\ 



(A.41) 



(A.42) 



(A.43) 



a. The hyperfine contribution 



The tensor amplitude of the hyperfine interaction is obtained by separating the spin dependence in eq. ( |A.35b . This means, 
one has to evaluate 



92 



d(Ji,idah,j 



SS,2 



l'=0 1=0 



X n 



[0] 



(p',]5). 



(A.44) 



It suffices to consider dfj'^ for the special case, for which the scattering plane coincides with the x-z-plane, i.e. p' 
(sin 6*, 0, cos 6). First one notes that then 



'{l-m)l ( V I 2 



[l + m)\ \m 



P,™ (cos 61), 



(A.45) 
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where denotes the associated Legendre function, and thus 
92 



(A.47) 



dcri^idcThj 

where 

With these expressions one finds for the tensor part of the hyperfine contribution 



where for to = 0, 1, 2 



Explicitly, with 



j[2],0 „ ^SS 

-'is ~ 



where I have introduced for convenience 



(A.48) 



fc=|i-2| \ / \ / 



'^11/22 - --y^ ± > 1^33 — '^"lOy 3' 0-12 — ^y'>|m|2 , Cri3 " " y ^|m|l . Cr23 — "2 1™!^' (.A.5U) 

one obtains for the nonvanishing components 

00 

421,0 ^ e«^^5PP,(cose), (A.51) 

00 ^ 

41/22 = ±^"" Y.S\^Pl(co.B) - i4f °, (A.52) 

00 

^5'/3P;i(cos6') , (A.53) 



C33 _ 1 00 oil _ 1 3a-2)! _ 1 m-\)\ 1 

-2*" -2Y2(7T2)! " ' ""2 V2(ITT)! ' ■ ^ ^ ^ 

It is useful to separate the yy-independent contributions, constituting the plane wave approximation. One finds explicitly the 
following detailed expressions 

41'°('?) = ( sin2(^/2) - 3 + E ^''(^) (COS ^)) ' (A.55) 

;=o 

^ 00 

41'°W = c«^(-2sinW + E^'"(^)^''('=°s^))' (^-57) 

i=i 

where the coefficients (77) vanish for ry = 0. In detail one finds for i = j = 3 and Z = 0, 1 



533/ ^) ^ _3 r l + i??/i(?7) _ i77(5 - Z7?) ^!!Lfi ^ 

^ 5I-I-W 3 6(2- w){3-«r)) 1-mV 



1 — if] 3 6(2 — i?7)(3 — «?7) 1 — iry V {2 — ir]){3 — irj) J . 



(A.59) 
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and for Z > 1 



(2l-l)(2l + 3) 



Here I have introduced for convenience 



biiv) = 



2l-l + ir] 



{21 l)l{l - 1 + ir]){l + i-q) ' 



One should note that the coefficients Sf^ behave as 1/P for ^ — > oo. 



For i — j ~ 1 one obtains (note I > 1) 

(2; + l)(l-e2**') 



2(Z-l)/(; + !)(/ + 2) 



(^(&Kr/)-6;+2(r/)*) 



2l{l + l){2l-l){2l + S) 



The coefficient S}^ behaves as / for I — > oo. 



Finally, for 5/ one obtains for / = 1 



3 
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and for Z > 1 



1(1 + l)i2l - 1){21 + 3) 



(A.60) 



(A.61) 



(A.62) 



(A.63) 



(A.64) 



The coefficient Sj^ behaves like / for / oo and vanishes for ry = 0. 

The convergence of the partial wave series is quite good in general as is demonstrated in Fig. lAll for r; = 2. Only dj^' '"(77, ( 
shows a slower convergence at very small angles. 



b. The spin-orbit contribution 

Acording to eq. ( IA.34l i the spin-orbit strength is given by 

00 



1=1 



For the chosen reference frame, i.e. p — (0,0,1) and p ' — (cos </> sin 0, sin sin 6*, cos Q), one obtains 



- sm( 
cose 



^viAfl.p) = 0, 
and thus the spin-orbit vector b has the form 



, P X P 

p' X p 



with 



-i 00 

&o('7,e) = --c^^^/3/(?7)Fi^cos( 



i=\ 



where 



1 « - e^*'^' 
AW = -;^e2-'i?„ = ^^^(2^ + 1 + /,(,?)). 



(A.65) 



(A.66) 
(A.67) 

(A.68) 
(A.69) 
(A.70) 
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This form is not well suited for a numerical evaluation, because even for 77 = the sum extends up to infinity. Therefore, it is 
more advantageous to separate the 77-independent part writing 



21 + 1 
l{l + l) 



with 



ft" 



{2l + l){e'^^' -l) + fi{r,) 



The coefficient /3,'' vanishes for 77 = 0. For the ij independent part one can rearrange the sum by using 



and one finds 



(A.71) 



(A.72) 



(A.73) 



1 



1=2 1=0 



Vl~x^ 



Then one obtains for bo{ri, 9) 



ic^«(cot(0/2)-5]/3,''P,i( 



(A.74) 



(A.75) 



yielding 60 (0, 6) in accordance with eq. ( l23T l. One can evaluate directly this expression or rearrange also the remaining sum 
yielding 



where for / = 0, 1 



1 1 °° 

- c^^ ( cot (^/2) - — 5] ei (77) Pi (cos ( 

1=0 



(A.76) 



(A.77) 



and for Z > 1 



{i + m + 2) (1-1)1 ^, 

2/ + 3 ^'+' 2l-l'^-' 

2^ + 1 fiiv) 



2 e^**' 



/ . ( 21-1 2l + -6]\ 



177 



21 + 3- 



{l + iTj)il + 1 - irj) {2l-l){2l + 3) 
1 1 



.{2l + 'S){l + l- iT])^ {21 - 1){1 + irjY 



(A.78) 



For I — > 00 the coefficients /3;' and e; behave as 1// resulting in a considerably slower convergence than for the hyperfine 



amplitude. This is demonstrated for the region of small angles and near 9 = 180° in Fig. I A2 1 for r/ = 2 for various ^ma^-values 
up to Imax = 5000 with the expansion of eq. ( |A.76l l. For comparison the result for the integral representation is also shown. The 
expansion of eq. ( IA.75b gives for large absolute values of the same result. However, for small 77 it results in small oscillations 
around the result of the other expansion at small angles. 



3. Comparison of the two methods 



For the hyperfine and spin-orbit interactions both methods give identical results for the corresponding amplitudes as is demon- 
strated for the hyperfine amplitudes in Figs. IA3 I through lAS I for -q — ±2. The above mentioned numerical problem which arises in 
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the integral representation for higher negative 77 is illustrated in Fig. lA6l where. as an example, the reduced amplitudes TZ§^ {rj, 9) 
and I^^{ri, 9) for 77 = 2.5 are plotted. The two curves labeled "IR" and "PWE" are almost indistinguishable. This means also 

complete agreement for the amplitudes for positive 7/ > as shown in the lower two panels of Fig. lATl forthe amplitude dgg {rj, 9). 
Furthermore, one notes that for this 77-value the two reduced amplitudes in Fig. I A6 1 are almost equal to a high degree of accuracy. 
Thus, this feature creates the numerical problem mentioned above for the integral representation method, because for negative 
values the amplitudes are represented as differences of the reduced amplitudes (see eqs. (lA.lSb and (IA.24b ). Consequently, large 
cancellations occur which reduce the numerical accuracy more and more with inreasing absolute values for negative 77. In fact, 
accuracy is lost for about —rj > 2. This is demonstrated in the upper two panels of Fig. |A7| for rj = —2.5 where one readily 
notes the onset of some numerical instabilities for the curves labeled "IR", in particular at small and large angles. This limits at 
present the numerical application of the integral representation method. 

A comparison of the two methods for the spin-orbit amplitude is shown in Fig. IA8I again for 77 = ±2, where one readily 
notes very good agreement. The same numerical problem of the integral representation method arises also in this case for large 
negative rj. 
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FIG. Al: Convergence of the partial wave expansion (PWE) of the hyperfine amplitudes 6) for various Imax as indicated in the legend 

for ?7 = 2 for (ij) = (11) (upper panels), (ij) = (13) (middle panels), and (ij) = (33) (lower panels). For comparison the result of the 
integral representation (IR, solid curves) is also shown. 
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FIG. A2: Convergence of the partial wave expansion (PWE) of the spin-orbit amplitude b''{rj, 9) for various Imax as indicated in the legend 
for 77 = 2 near 6 = 0° (upper panels) and near 180° (lower panels). For comparison the result of the integral representation (IR, solid curves) 
is also shown. 
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FIG. A6: Reduced amplitudes TZH {ri, 0) and I33 (r?, 0) for r? = 2.5 for the integral representation (IR) and the partial wave expansion (PWE). 
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FIG. A7: Hyperfine amplitude d33''^('?, 0) for rj = ±2.5 for the integral representation (IR) and the partial wave expansion (PWE). 
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FIG. A8: Spin-orbit amplitude 6o(??, 0) for r? = ±2 for the integral representation (IR) and the partial wave expansion (PWE). 



